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Abstract 

It is shown that irreversible classical cellular automata can be performed 
by quantum algorithm using additional ancilla registers. The algorithm for 
cellular automata states analysis has been proposed. This algorithm is based 
on the elements of Grover's algorithm - the inversion of amplitude of searched 
states and unitary transform of inversion about the average. The inversion of 
searched states amplitudes can be performed by quantum Toffoli gate. 

Introduction 

The cellular automata (CA) is a computational model which has been studied in 
many scientific works . These models are based on simple updating rules and show 
the perspective for computational applications. 

The aim of our study is to design the algorithms which allow to use quantum 
parallelism for the investigation of CA evolution for all initial states simultaneously. 
For example, if one dimensional CA consists from only 100 cells with 2 possible 
states, then one quantum evolution can evaluate 2^°° evolutions for the same number 
of initial states. It is important to answer a question - does or does not some searched 
finish state ("good state") in CA evolution for analyzed CA updating rules exist? We 
will use Grover's algorithms ideas to answer this question in case of one dimensional 
quantum cellular automata. 

Quantum cellular automata (QCA) are being investigated in many modern works. 
In the work [2] QCA for universal quantum computation has been presented. In [2111] 
one dimensional QCA has been investigated. In [S] QCA formalism based on a lat- 
tice of qubits has been described. In [6, 2] reversible QCA was considered. In [8l |9j 
the survey of QCA investigations has been presented. 



One Dimension Cellular Automata 

Let us consider the ability of implementation of synchronous cellular automata by 
the quantum logical gates when all cells switch to a new state simultaneously. Such a 
circuit can be performed by using additinal qubits named ancilla. Let us consider for 
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determinance one of CA rules for updating CA states. For example, a cell switches 
to state "0" if its neighbors have the same value, otherwise this cell switches to value 
" 1" . This rule can be written as 

x'[i\ ■= x[i - 1] © x[i + 1]. (1) 

In the classical CA such evolution is irreversible. If the results of previous iteration 
are saved in an additional ancilla, then CA evolution reversibility and unitarity can 
be achieved. Consider unitary operator for performing rule ([T]). This operator acts 
on three qubits - |xi_i), , \ai) two neighbors and ancilla from the additional 

register: 

C = {X ® X ® I) ■ T ■ {X ^ X ^ I) ■ T ■ {I ® I ® X). 
Before the iteration all additional ancilla are in the state 



|0,0,...,0)„ 

The series of unitary transformations for qubits can be written as 



Ci 
Ci 

a 



\Q)i-lAQ)^+lAQ)a,t+l -> |0)i_l,t|0),+i,t|0),,t+i 
|l)i-l,t|l)i+l,i|0)a,t+l |l)j-l,t|l)i+l,t|0)a,t+l 

a,t+l 

|0)i_i,t|l)i+i,t|0)a,t+l \Q)^-lAl)^+l,t\l)a,t+l 



(2) 
(3) 



d: I ® I ® X\xi-i, Zi+i, 0) Xi+i, 1) \xi^i,Xi+i, 1 © Xi^iXi^i) 

^ X © a: © Xj+i, 1 © Xi^iXi+i) hxi_i, -^Xi+i, 1 © Xi^iXi+i) 

— )■ T|— "iXi^i, 1 © Xi^iXi^i) — )• — iXj-i-i, 1 © Xi^iXi^i © 

A © A © I\^Xi^i, -^Xi+i, 1 © Xi^iXi+i © -.Xj-i-iXi+i) \xi-i,Xi+i, (4) 

The action of operator Cj on the possible qubits values can be described as 



(5) 



CA updating rules dlJ-dS]) can be described by quantum circuit in figJU 
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Figure 1: Quantum circuit implementation for CA updating rules 
CA state at j iteration can be described using the following operator 

A, = C^...C2C^. (6) 
Let us consider the first CA iteration with the additional ancilla register. 

A, : |x)o|0)i ^ \x)o\f{x))i. (7) 
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where |x) = | ) - is initial CA register, f{x) - is some function which 

expresses CA rules for CA updating ([I])-®. On the next iteration the new ancilla 
register is being added. It can be described as 

A2 : \x)o\f{x))i\0)T -> |x)o|/(x))i|/(/(x)))2. (8) 
CA evolution after m iteration can be written as 

A^... A,A, : |x)o|0)r . . . |0)r ^ \x)o\fix)), . . . If^^^Hx))^. (9) 

Our next step is to perform the iteration which is inversed to the one that is defined 
by operator Am~i, where 

Am-Mm-l)-' = 1. (10) 

Applying of such an operator to the registers system after m iterations ([9]) can be 
written as 

(A^.,)-' : |x)o|/(x))i . . . \f^'^-'\x))^^^\f'^\x))^ ^ 

^\x)o\f{x)h...\0)T-i\f^"'\x)U. (11) 

A set of inversed transformations can be written as 

{A^)-' . . . (A„_2)-^(A„-i)-^ : \x)o\f{x))^ . . . \f^-'\x))m-i\f^"'\x))n. -> 

->l^)o|0)r---|0)ril/^'"n^))^- (12) 

As a result of applying the operators of inversed evolutions, the registers of additional 
qubits will be in initial states and can be removed without the affection on other 
qubits. General evolution of CA can be described by operator 

UcA = (Ai)-^ . . . {Am-2)'\A„,^i)-^A^A^.i . . . Ai. (13) 

The resulted evolution of quantum CA can be written as 

f/cA:|x)|0)®"^|x)|/('")(x)). (14) 

Let us consider the steps of implementation of quantum cellular automata: 

1. At the first step the CA register in the basic state is initialized 

\XI,X2, ■ ■ ■ Xn)n |0l, O2, • • • 0„)„ (15) 

2. Let us apply the Hadamar operator to each qubit of CA register. This operator 
is given by matrix 

-^(;-\)- 

3. As a result we obtain the following superposition 

1 x=2" 

i/^'^|0i,02,...0„)„ = -= k) (17) 
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where |a;) denotes computational basis states \xi,X2, ■ ■ ■Xn)n 

The Hadamar transformation generates the superposition of all possible 2" 
states with the equal amplitudes. 

Let us apply unitary transformation ([H]) 

\^ca) = UcA (i/^"|0)®") 10)®'^ = 4^ i: (18) 

Superposition I^E'ca) includes all 2" results of evolutions initial states of CA 
with dimension n. With one quantum evolution of CA, all possible CA evolu- 
tions for initial states with defined CA updating rules are performed. 



Analysis of Cellular Automata Evolution 

Let us consider the abihty to amplify the amplitudes of searched states with using 
the ideas of Grover's algorithm which is being being used for searching files in the 
quantum database pUl ITT| IT^ . The difference between this task and Grover's al- 
gorithm is that unknown states are not being searched in this task, these searched 
states are known and we need to get the answer that such states exist in CA evo- 
lution. Let us introduce an additional qubit - ancilla which will be controlled by 
n-qubit Toffoli gate where n-qubit of CA state are controlling the ancilla state in 
the Toffoli gate. This gate can be described by unitary matrix 
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(19) 



Let us consider the unitary operator which is defined as tensor product of one 
qubits operators 



i=l 



where 
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(20) 



(21) 



The operator St flips searched states |ei, 62, . . . e„) into the state |li, I2, • • • In)- It 
is necessary for implementation of anclilla inversion for searched states by Toffoli 
transformation f lTOj) Apply hadamar operator to ancilla \z) in prepared state |1) 



\z) = H\l) 



V2 



(|o)-|i)) 



(22) 



Let consider an operator 



Ut = (/^" ® 5r ® /) (/^" ® T„) (/«" ^St®I) (23) 

Let apply this operator to the system of qubits registers \x)\f^"^\x))\z). the first 
group of operators on the right side which is delimited by brackets performs flipping 
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into the state |li, I2, . . . In) of searched finish states, the second group performs 
ancilla \z) inversion for searched states and the third group returns states which are 
changed by the first group into the state before applying the operator. Additional 
controlled qubit stays in the state fl22|l . The result of the operator fl23l) acting is 



[/t(|x)|/(™)(x))|^)) = ^ K: e (24) 

where Xg - is set of initial states \x), which leads to the searched states \q) in the CA 
evolutional process after m CA iterations. Ancilla I2:) in a new basis leaves in the 
unchanged state, but in the superpositions states the inversion of amplitudes signs 
for subsystem J2xeXq \^)\q) appears. It is caused by the transformation of ancilla \z) 
to new state (1221) before applying the operator (^^. 

Let us consider the inversion operator from Grover's algorithm 

t/G = 2|^e)(*c|-/ (25) 

where 

1 1=2" 

ivi/,) = Hnor"" = ^ E N) (26) 

V^" i=i 

This operator will be used for the states J2xeXq 1^)11) amplitudes amplification. The 
operator Ug refiect any vector around the axis defined by vector |\E'c)- The operator 
Ug can be decomposed by one-qubit operators 

Ug = /f®'^(2|0)(0| - J)//®" (27) 

The operator Ug also has been called as the operator of inversion about the average. 
The transformation (1251) can also be defined as 

f/c : E«iK) ^E(2^-ai)K) (28) 

i i 

where A - means the average value of Oj amplitudes. 

If searched state Z^zex, \^)\q) can be realized only once from one initial CA state, 
then its amplitude is 

/3 = ^ (29) 



It can be shown that performing of CA iteration 

UGUT{\x)\f^-\x))\z)) (30) 

we obtain the amplitude amplification by 3 times that match the implementation of 
inversion iteration in the Grover's algorithm (TUl [HI [12] . The transformation UgUt 
can be described by quantum circuit in figJ2J 

Let us consider the number of iterations necessary for enough amplification of 
searched states amplitudes. If only one initial state \x) of cellular automata leads 
to the searched finish state then using the ideas similar to Grover's algorithms 
[TUl [TT| [T2] we can find the optimal number of unitary transform UgUt ■ 

k^^VN,N = 2^^ (31) 
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Figure 2: Quantum circuit implementation of the transformation UqUt 



It follows from this result that algorithm complexity is 0(v A^). In comparision with 
analogical classical algorithm 0{N) it means polinomial computational speedup. If 
instead 1 matching initial state there are / matching states which lead to the searched 
finish state, then 

' ^ (32) 



k 



4V / 



if / is unknown, then Grover's algorithm could be run several times when 



L = l,2,4,: 



(33) 



It can be shown, that for such series calculation the complexity will be still 0{\/N). 



Conclusion 

It is shown that irreversible classical cellular automata can be performed by quan- 
tum algorithm using additional ancilla registers. The algorithm for cellular au- 
tomata states analysis has been proposed. This algorithm is based on the elements 
of Grover's algorithm - the inversion of amplitude of searched states and unitary 
transform of inversion about the average. The inversion of searched states ampli- 
tudes can be performed by quantum Toffoli gate. 
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